Abstract. A sample of OB stars (eleven Of, one O and 
one B supergiant) has been surveyed with the Australia 
Telescope Compact Array at 4.8 and 8.64 GHz with a 
resolution of ~2"-4". Five stars were detected; three of 
them have negative spectral indices, consistent with non- 
thermal emission, and two have positive indices. The ther- 
mal radiation from HD 150135 and HD 163181 can be 
explained as coming from an optically thick ionized stel- 
lar wind. The non-thermal radiation from CD-47° 4551, 
' J 1 ■ HD 124314 and HD 150136 possibly comes from strong 
[ shocks in the wind itself and/or in the wind colliding re- 

■ gion if the stars have a massive early- type companion. 
' The percentage of non-thermal emitters among detected 

j>^> O stars has increased up to ~50%. The Of star HD 124314 
d ' clearly shows flux density variations. Mass loss rates (or 

■ upper limits) were derived for all the observed stars and 
[ the results compared with non-radio measurements and 

■ theoretical predictions. 

' Key virords: stars: early-type - stars: mass loss - stars: 
^sD winds, outflows - radio continuum: stars 
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1. Introduction 

Of stars are characterized by N[iii] A A 4630-34 A emission 
Hnes and He [11] A 4686 A emission or absorption lines in 
their spectra. According to Maeder (1990) a massive O 
star evolves into an Of star which, in turn, goes through 
the blue supergiant and luminous blue variable phases, 
and after that becomes a Wolf Rayet (WR) star. In ad- 
dition, Of stars have very strong winds with gas termi- 
nal velocities up to 3000 km s~^ (Chlebowski & Garmany 
1991). The velocities, derived from P Cygni profiles, are 
appreciably larger than the stellar escape velocity, prov- 
ing that important mass loss is taking place in Of stars. 
Typical values for the mass loss rate, M, of Of stars are 
<, IQ-^ Mgyr-i (Lamers et al. 1995). 

Like for WR stars, the action of the wind on the am- 
bient interstellar medium produces a cavity of hot rar- 
efied gas surrounded by a slowly expanding envelope. The 
shells can be detected at optical and radio wavelengths 
(Lozinskaya 1982, Cappa & Benaglia 1998, Benaglia & 
Cappa 1999). Using high resolution, the wind region itself 
is observable in the radio continuum emission. The stellar 
mass loss rate (or an upper limit) can be estimated from 
the radio flux density. Accurate determinations of M are 
essential in modeling the stellar evolution, describing the 
Of-stellar wind and changes in the ambient gas caused by 
the mass outflow. Since radio-determined mass loss rates 
turned out to be among the most reliable ones, they are 
used to calibrate M obtained with other methods, i.e. in- 
volving Ha equivalent widths or UV profiles. 

It is commonly assumed that the dominant driving 
mechanism for the winds is radiation pressure. However, 
M values predicted by radiation-driven wind models differ 
up to a factor of two with the observed mass loss rates for 
0-type supergiants (Vink et al. 2000). Puis et al. (1996) 
found the discrepancies between observed and predicted 
mass loss rates diminished when considering more accu- 
rate input parameters of the involved atomic physics (see 
Kudritzki & Puis 2000 for a comprehensive review). In 
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this respect, new results on M are crucial for testing the 
existing wind models and theoretical predictions. 

Using the Very Large Array (VLA), systematic obser- 
vations of O, B and WR stars were carried out by Abbott 
et al. (1986) and Bieging et al. (1989), while more re- 
cently Scuderi et al. (1998) also observed Of as well as O 
and B-type stars. With the Australia Telescope Compact 
Array (ATCAQ), Leitherer, Chapman & Koribalski (1997) 
extended the study towards southern declinations, observ- 
ing a complete sample of WR stars closer than 3 kpc. Up 
to now, only a few southern OB stars were observed (e.g. 
Leitherer, Chapman & Koribalski 1995). We have initiated 
a programme with the aim of detecting bright early-type 
stars that remain unobserved at radio wavelengths in the 
southern sky. Here we present the results obtained so far, 
after the first two observing campaigns. 

The next section gives characteristics of the selected 
sample of OB stars. Section 3 describes the observations 
and reduction procedure. The theoretical formulae used 
here are briefiy described in Section 4. The results are 
presented in Section 5 and discussed in Section 6. We close 
with the conclusions in Section 7. 

2. The target stars 

From the lists published by Cruz-Gonzalez et al. (1974) 
and Garmany et al. (1982) we selected all southern stars 
(declinations south of -24°), closer than 3 kpc and still un- 
observed at cm wavelengths. These stars form our working 
database. 

We derived an expected mass loss rate with the approx- 
imation given by Lamers & Leitherer (1993). Their for- 
mula gives M as a function of effective temperature, T^s-, 
and stellar luminosity, i». The values for Tcff and L, were 
taken from Vacca et al. (1996). By means of the standard 
model for thermal emission from stellar winds (Wright & 
Barlow 1975, Panagia & Felli 1975) we computed a min- 
imum detectable mass loss rate, assuming the radio flux 
as two times the predicted r.m.s. after 3 hours integration 

^ The Australia Telescope is funded by the Commonwealth 
of Australia for operation as a National Facility managed by 
CSIRO. 
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time with the ATCA. The predicted r.ni.s. noise resulted 
in 0.050 mjy at 3 cm, and 0.054 mJy at 6 cm, if contin- 
uum emission from the target stars was not contaminated 
with confusion from other sources in the field. Finally, we 
defined an index F equal to the expected mass loss rate 
over the minimum detectable mass loss rate, and ranked 
our stars according to that index. 

The present target stars were chosen from the database 
depending on the allocated time, and favoring the objects 
with highest index F, i.e. the closest and earliest. The 
candidates are listed in Table 1. For this sample, the min- 
imum detectable mass loss rates were (1 - 3) 10~^ 
yr-i. 

We have included the already observed star HD 57060 
(29 CMa), which was not detected with the VLA (Abbott 
et al. 1980), in the target list because, besides matching 
most of the criteria, it is a very interesting Of-|-0 binary 
system. The r.m.s. noise expected from the ATCA, 0.050 
mJy, is 8(7 lower than the reported VLA flux limit of 0.4 
mJy. The masses of both the primary and secondary stars 
are in the range 20-30 M©, the primary being the more 
evolved star (Hutchings 1977, Stickland 1989). Wiggs & 
Gies (1993) interpreted the Ha P Cygni profile by means 
of a model where emission comes from inside the wind of 
the primary and from a plane midway between the stars 
where the winds collide. 

Ha emission has also been reported in HD 163181 
(Thaller 1997) and thus this binary star can be suspected 
to have colliding winds. 

VLA 6-cm observations towards the field of HDE 
319718 (N35 according to Neckel 1984) were presented by 
Felh et al. (1990). Again, using the ATCA at 3 and 6 cm 
we intended to obtain a more sensitive image. 

Table 1 lists the stellar parameters, beginning with the 
name, spectral classification and photometric data. Infor- 
mation concerning possible membership to a stellar as- 
sociation is given in column 6. The spectro-photomctric 
distance was estimated using the and [B — V)o values 
from Vacca et al. (1996) and Schmidt-Kaler (1982), re- 
spectively. In the case of HD 163181, the absolute magni- 
tude was extrapolated. In all cases but HD 101205, the dis- 
tances estimated by us were, within errors, in agreement 
with the distances to the OB associations. We adopted 
as the stellar distances d the cluster distances whenever 
available, and the spectro-photometric distances for field 
stars. The mass loss rates quoted in Table 1 were derived 
by means of optical and UV spectra. When no measured 
terminal velocities [voo) were found in the literature, the 
values were interpolated from Prinja et al. (1990) and ap- 
pear in brackets. Column 9 accounts for the binary status. 
The label "BIN" stands for the stars catalogued as bina- 
ries with computed orbits. Following Gies (1987) SB2? 
means that double lines have been reported, but no orbit 
was determined; SBl? indicates that the star is a possible 
spectroscopic binary; and C, that the radial velocity was 
taken as constant, i.e. as for single stars. 



3. Observations and analysis 

The observations presented here were obtained with the 
Australia Telescope Compact Array in two campaigns: 1) 
in 1998, February 16/17, using the 6 A configuration, and 
2) in 2000, March 30 + April 1/2, with the 6D configu- 
ration, collecting a total of 52 h observing time. In both 
sessions the maximum baseline was 6 km. 

Each source was observed simultaneously at two fre- 
quencies: 8.64 GHz (3 cm) and 4.8 GHz (6 cm), with a 
total bandwidth each of 128 MHz over 32 channels. The 
average time on source was 2.5 to 3 h, exception made 
for HD 124314 which was observed on three occasions, in 
looking for flux variability. The stars were observed dur- 
ing intervals of 15 or 30 minutes, interleaved with short 
observations of close phase calibrators, and spanned over 
an LST range of ~11 hours for most of them. The allo- 
cated time allowed us to track the sources HD 57060, HD 
163181 and HDE 319718 for -7 h each. The flux density 
scale was calibrated using observations of the primary cal- 
ibrator, 1934-638, assuming flux densities of 2.84 Jy at 3 
cm and 5.83 Jy at 6 cm. 

The data were reduced with Miriad routines and an- 
alyzed with the AIPS package. After data calibration the 
visibilities were Fourier-transformed using 'natural-' and 
'uniform'-weighting. The images obtained in the former 
way turned out to have better signal to noise ratio. The 
restoring beams were J^l'.'5 at 3 cm and ^3" at 6 cm. In 
many cases, the diffuse emission from extended sources 
towards the observed fields had to be removed by taking 
out the visibilities corresponding to the shortest baselines. 

Flux density values of the thirteen observed sources are 
listed in Tables 2a and 2b. The three entries of HD 124314 
correspond to the different observing dates (CI: Campaign 
1, and C2: Campaign 2). The optical positions were taken 
from the SIMBAD database. Four sources were detected 
at both frequencies, with flux densities greater than 3cr. 
HD 150135 was detected at 3 cm only; at 6 cm it appears 
at a 2.5cr level. All detected sources were unresolved: CD- 
47° 4551, HD 124314, HD 150135, HD 150136 and HD 
163181. Bi-dimensional Gaussians were fitted to the five 
sources using Miriad task ImRt, and the positions of the 
maxima are given as the radio position in Table 2a. The 
differences with optical coordinates, whenever present, are 
< 0'.'5, so the optical and radio positions are in good agree- 
ment within the errors. The flux density was computed 
after correcting the zero level in the images. The flux den- 
sity error in the Tables is equal to the image r.m.s. noise 
(cr). For the undetected stars, a maximum flux density of 
3(7 is quoted. The integration time on source is listed in 
Tables 2a and 2b as tint. 

Figures 1 to 4 display the detected stars, at 3 and 6 
cm. The restoring beam is plotted in each case. 



p. Benaglia et al.: Mass loss rate of OB stars 



o 



—I 



-484538 

40 
42 
44 
46 
48 
50 
52 
54 
56 



a. HD150136 and HD150135, 3cm 





o 




164121.0 20.5 20.0 19.5 19.0 164121.0 
RIGHT ASCENSION (J2000) 



20.5 20.0 19.5 19.0 
RIGHT ASCENSION (J2000) 



Fig. 3. a) Contour map of HD 150135 and HD 150136 at 3 cm. Contour levels are -0.08, 0.12(4cr), 0.2, 0.3, 0.6, 1.0, 
1.4, 1.8 and 2.2 mJy per beam area. The restoring beam, plotted at the bottom left corner, is 2.6x1.4 arcsec in P.A. 
-33°. b) The same as a) at 6 cm. Contour levels are -0.15, 0.15(5c7), 0.2, 0.4, 1.0, 1.75, 2.5, 3.25, 4.0 and 4.75 mJy per 
beam area. The restoring beam is 4.1x2.0 arcsec in P.A. -33°. 



4. Radio emission from stelleir winds 

Stars with ionized winds show a flux density excess to- 
wards long wavelengths. Assuming the continuum radia- 
tion is due to free- free emission (thermal brcmsstrahlung) , 
Wright & Barlow (1975) and Panagia & FeUi (1975) de- 
duced the spectrum in the radio and IR energy ranges. 
Their model was developed for a uniform flow expanding 
at constant velocity, and allows to derive the stellar mass 
loss rate as a function of the measured flux density as: 

c3/4 ,3/2 

M = 5.32xl0-^^V^-^^ MQyr-\ (1) 

Here d is the stellar distance in kpc, Woo the wind ter- 
minal velocity in km s~^, the flux density in mJy; fi 
stands for the mean molecular weight, 7 is the mean num- 
ber of electrons per ion and Z, the rms ionic charge. 

In the range of the wind electron temperatures Te and 
observing frequencies we are dealing with, the free-free 
Gaunt factor can be approximated by: 

/ 2^3/2 \ 

5^ = 9.77 I l + 0.131og^J . 

The model predicts a spectral index a = 0.6 (5^ ex :/") 
in accordance with measurements of various WR and O 
stars. If the spectral index differs from 0.6, Eq. (1) yields 
an upper limit to the stellar mass loss rate. However, when 
detected at two or more radio continuum frequencies, an 
important number of stars show spectral indices approach- 
ing zero and even negative values. The interpretation of 
these results is consistent with a picture where the ob- 
served flux density has an important non-thermal contri- 
bution. The main non-thermal process which is relevant 



for early type stars is synchrotron radiation by relativis- 
tic electrons, accelerated by first order Fermi mechanism 
at strong shock fronts. The presence of non-thermal ra- 
dio emission provides indirect evidence for the presence 
of magnetic fields. In some binary stellar systems the 
electrons seem to be accelerated at the wind collision re- 
gion (i.e. Cyg 0B2 No. 5 and WR 140, Contreras et al. 
1997, 1996). Eichler & Usov (1993) have shown that stel- 
lar winds collision in early type binaries may be strong 
sources of non-thermal radio and 7-ray emission (see also 
Benaglia et al. 2001 for an updated discussion). 

For single stars, White (1985) has shown that electrons 
can be accelerated up to relativistic energies at the mul- 
tiple shocks formed at the base of the wind by line-driven 
instabilities. The derived radio spectrum follows an index 
of ~0.5 up to a turnover around a few GHz. The relativis- 
tic electron population would generate inverse-Compton 
7-rays too. Benaglia et al. (2001) have estimated the ex- 
pected 7-ray luminosity from the wind colliding region, 
the terminal shock of the wind and the zone of the base 
of the wind, for Cyg 0B2 No. 5, using the radio results of 
Contreras et al. (1997). 

An interesting feature seen in some non-thermal 
sources is their variations in radio luminosity with time. 
Abbott et al. (1984) found this behaviour in 9Sgr and 
Cyg 0B2 No. 9. Cyg 0B2 No. 5 has been extensively mon- 
itored, since ~1980 (see Contreras et al. 1997 and refer- 
ences therein). The radio source seems to switch between 
states of high and low emission with a period of 7 years. 
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Fig. 1. a) Contour map of CD-47° 4551 at 3 cm. Contour 

levels arc -0.1, 0.12(2.5cr), 0.18, 0.3, 0.5, 0.7, 0.95 and 1.2 
mJy per beam area. The restoring beam, plotted at the 
bottom left corner, is 2.7x1.0 arcsec in P.A. -2°. b) The 
same as a), at 6 cm. Contour levels arc 0.1, 0.18(3.5(t), 
0.3, 0.5, 0.7, 1.2, 1.7 and 2.2 mJy per beam area. The 
restoring beam is 4.5x1.6 arcsec in P.A. -2°. Negative con- 
tours are dashed. The optical position of the star is marked 
by a cross. 



5. Results 

Wc detected five stars out of thirteen: CD-47° 4551, HD 
124314, HD 150135, HD 150136 and HD 163181. Table 3 
presents the spectral indices. The radio emission of the 
candidates CD-47°4551, HD 124314 and HD 150136 is 
characterized by a negative spectral index between 4800 
and 8640 MHz. The one of HD 124314 displays strong 
variations. The flux densities of these stars seem arising 
mainly of a non-thermal process. The only B star in the 
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Fig. 2. a) Contour map of HD 124314 at 3 cm. Contour 
levels are -0.15, 0.15(3(7), 0.5, 1.0, 1.5, 2.0 and 2.5 mJy 

per beam area. The restoring beam, plotted at the upper 
right corner, is 1.9x1.2 arcsec in P.A. -52°. b) The same 
as a) at 6 cm. Contour levels are -0.15, 0.15(3(t), 0.6, 1.2, 
1.7, 2.2, 2.7 and 3.2 mJy per beam area. The restoring 
beam is 3.2x2.1 arcsec in P.A. -52°. 



group, HD 163181, and HD 150135 show positive spectral 
indices, possibly indicating thermal emission as dominant. 

From hereafter we shall call "thermal" sources the ones 
with positive spectral index, and "non-thermal" if the in- 
dex is negative. 

The measured flux densities at 3 cm were used to de- 
rive the mass loss rates of all target stars, because of their 
smaller errors compared to the flux densities at 6 cm. Ac- 
cording to their evolutionary status, for most of the target 
0-type stars the computed mean molecular weight was 
1.5. We adopted a value of 1.3 for HD 124314 and HD 
150135, the less evolved stars, and 1.6 for the only B-type 
star in the group. In deriving the rms ionic charge and 
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Table 3. Spectral indices of the detected stars 



1756 16.8 16.4 16.0 15.6 
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Fig. 4. a) Contour map of HD 163181 at 3 cm. Contour 
levels are -0.15, 0.15(3cr), 0.2, 0.25, 0.3, 0.35 and 0.4 mJy 
per beam area. The restoring beam, plotted at the upper 
right corner, is 3.9x1.1 arcsec in P.A. -34°. b) The same 
as a) at 6 cm. Contour levels are -0.13, 0.13(2.5cr), 0.17, 
0.21 and 0.25 mJy per beam area. The restoring beam is 
6.4x1.7 arcsec in P.A. -34°. 



the mean number of electrons per ion wc assiimcd that 
both Hydrogen and Helium were singly ionized, and that 
in the region where the radio emission comes from, the 
He+"'" had the chance to recombine. The electron temper- 
ature was assumed as 0.4 the effective stellar temperature. 
The computed values of n, Z, and 7 are listed in Table 4, 
together with the free-free Gaunt factor at 3 cm. 

We derived the mass loss rates by means of Eq. (1), 
independently of the behaviour of the flux density with fre- 
quency (see Table 4) . The results for candidates with spec- 
tral index a < appear in parenthesis. For non-detections 



Star 


Spectral index 




Q3-6cm 


CD-47° 4551 


-0.89 ±0.06 


HD 124314 




CI: 02/16/1998 


-0.60 ± 0.07 


C2: 03/30/2000 


-0.87 ±0.06 


C2: 04/02/2000 


-0.46 ± 0.04 


HD 150135 


-hl.75±0.70 


HD 150136 


-1.29 ±0.03 


HD 163181 


+1.34 ±0.47 



and non- thermal sources, the derived mass loss rates must 
be considered as upper limits. 

In the computation of the mass loss rate errors we 
have adopted the relative deviations for /z, 7, gi, and Z of 
^10%, and for the terminal velocity an accuracy better 
than 20%. The flux density errors contribute with ~10% 
in the worst cases. Undoubtedly the main source of error 
in M is the uncertainty in stellar distances. If the relative 
error is of 20%, the resulting a-(log(M)) « 0.2. This last 
value increases to 0.3 in the more conservative case of 
distances accurate up to 30%. 

Inhomogeneities or anisotropy in the winds like 
dumpiness or coUimated flows are other causes of errors 
in M when estimated from the model of Wright & Barlow 
(1975). 

Among the eight non-detected stars, five had predicted 
mass loss rates of the order of the minimun detectable 
mass loss rates, and we expected that they would be de- 
tected at a 2(7 level. The stars HD 57060 and HD 112244 
have both expected mass loss rates, and mass loss rates 
derived by optical means appreciably larger than the min- 
imum detectable mass loss rates. 



5.1. Comments on individual stars 

HD 57060. This source was observed for 3 hours with 

the ATCA and was not detected. The r.m.s. of the final 
images agrees with the expected r.m.s. of 0.05 mJy at 3 
cm, and 0.054 mJy at 6 cm. According to these values we 
conclude that the flux density of this star must be below 
0.15 mJy, and the mass loss rate less than 3 x 10~^ Mq 
yr~^. Abbott et al. (1980) reported no detection at 6 cm 
above 0.4 mJy after 165 minutes of integration time with 
the VLA. 

CD- 4r 4551. Besides the target star, at 4800 MHz a 
second source appears 5.4" to the west, with a maximum 
over 6a- (Fig. lb). Its flux density is 0.33 mJy at 6 cm. 
The non detection of this small source at 3 cm means 
that the flux density is S^^m. < 0.15 mJy, letting us to 
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conclude that the spectral index is negative. At 1.7 kpc 
the spatial distance between the two maxima is about 
9000 AU. 

HD 94963, HD 97253, HD 101205, HD 112244, HD 
135240 and HD 135591. The observations allowed us 
to estimate upper limits to mass loss rates for these 
undetected stars. 

HD 124314. This star was first detected on February 
1998. The measured flux densities at 3 and 6 cm yield a 
non-thermal spectral index of —0.60. Fig. 2 depicts the 
stellar wind as observed the first day of C2. In order to 
look for variability with timescales of years and days, we 
reobserved the source on March 30 and April 2, 2000. 

For this particular star we produced, at each frequency, 
three images, maintaining the number of iterations, the 
baselines (greater than 10 kA) and the restoring beams 
(2" at 3 cm and 3.6" at 6 cm). Although we were not 
able to appreciate differences in shape between the maps 
made using data taken at different dates, the observed flux 
densities showed large variations (see Table 2a). 

At 3 cm the variability is appreciable between obser- 
vations of 1998 and 2000, as an increase of ~50%. Flux 
densities of C2 did not present differences within the error. 

The source displayed always variations when observed 
at 6 cm. The flux densities measured at 2000 are 65% 
and 40% greater than the corresponding one of 1998. Dur- 
ing the three days-interval of C2 the flux decreased about 
15%. 

The mass loss rate of HD 124314 was derived from 
the observations with the highest spectral index, i.e., for 
which the non-thermal contribution to the flux densities 
were the smallest. 

HD 150135 and HD 150136. Fig. 3 displays the radio 
images of these two stars. HD 150135 is barely detected 
at 6 cm, with Seem ~ 2.5cr. The earliest and strongest 
source, HD 150136, has a singular shape. Although 
elongated as the observing beam, it shows extensions to 
both sides, visible at the two observing frequencies. 

HD 163181. At 6 cm it is a 4ct detection, being the 
flux density about half of that at 3 cm (see Fig. 4 and 
Table 2b). 

HDE 319718. This star is located near the Hii region 
G353.2-h0.9. The 6-cm VLA map of the ionized gas (FcUi 
et al. 1990) has a peak of ^ 120 mJy, and covers a region 
of 2' in diameter, 1' north of our target star. To eliminate 
completely the emission from the extended sources, the 
ATCA image was built considering visibilities of baselines 
only with antenna 6. The candidate was undetected above 
0.15 mJy at 3 cm. However, upper limits to the flux density 
permitted us to derive a corresponding limit to the mass 
loss rate. 



5.2. Other sources in the fields 

Towards the direction of the target stars two previously 
unknown sources were detected, named SI and S2. They 
are visible at the two observed frequencies and are point 
sources for ATCA. 

SI is located in the field of HD 94963, at the ra- 
dio position a = 10''56'"29.27^ S ^ -61°43'1.19". The 
measured flux densities are S^cm = 3.8 ± 0.05 mJy and 
•^ecm = 6.4 ± 0.06 mJy, thus the spectral index results in 
-1.16. 

The second source is located at a = 15'^18'"35.49^ 
S = -60°29'59.44" in the field of HD 135591. Its radio 
flux densities are 530111 = 0.24 ± 0.05 mJy and Seem — 
0.60 ± 0.05 mJy, giving a spectral index of -0.68. 

Possible counterparts to these objects were searched 
using NED and SIMBAD databases and SkyView^ images 
from the Digitized Sky Surve yfl (DSS), with no positive 
results. The radio negative spectral indices indicate a non- 
thermal origin of the emission. This fact, together with the 
small angular size of the sources and the non existence of 
other related object suggest an extragalactic origin for SI 
and S2. 

6. Discussion 

6.1. Statistics on OB radio observations 

In order to perform some statistics on radio detections of 
OB stars, we have consulted the articles of Abbott et al. 
(1980, 1981, 1984), Bieging et al. (1989), Leitherer et al. 
(1995), Contreras et al. (1996), Waldron et al. (1998) and 
Scuderi et al. (1998), and added the results of this work. 
Regarding these investigations, about 120 OB stars have 
been observed, and less than 40 detected at least at one 
frequency. 

Excluding our observations, about 25 stars were de- 
tected at two or more frequencies, allowing the deter- 
mination of spectral indices. Of the almost 40 detected 
sources, ~50% classify as thermal, ~25% as non-thermal, 
and ^--^25% as composite or undecided. 

No non-thermal B stars were found, giving support to 
the suggestion of Bieging et al. (1989) that the synchrotron 
mechanism is less efficient at lower luminosities. 

After the detections presented here the thermal candi- 
dates stay around 50% between the detected O and B-type 
stars, while the non-thermal ones have increased from 25% 
to 30%. For detected 0-type stars, the figures are more 

^ SkyView was developed under NASA ADP Grant NAS5- 
32068 with Principal Investigator Thomas A. McGlynn under 
the auspices of the High Energy Astrophysics Science Archive 
Research Center (HEASARC) at the GSFC Laboratory for 
High Energy Astrophysics. 

The Digitized Sky Survey was produced at the Space Tele- 
scope Science Institute under US Government grant NAG W- 
2166. 
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striking: the thermal candidates are ^40% (before and af- 
ter this work), but the percentage of non-thermal O stars 
changes from 40% to almost 50% after the present results. 
Chapman et al. (1999) arrived at a similar result for WR 
stars, in surveying all candidates at negative declinations 
and distances less than 3 kpc. We think that the number 
of non-thermal candidates must be taken as a lower limit 
because: 

i) a single O star with powerful winds can radiate non- 
thermal radio emission due to relativistic leptons. Accord- 
ing to the model of White (1985), a spectral index of +0.5 
will be measured at low frequencies, mimicking a thermal 
spectrum. 

a) the flux density of a non-thermal source is known to 

vary, and at the time of the observation a minimum can 
be occurring, approaching a thermal spectrum. 

6.2. Flux variability and HD 124314 

It is important to analyze the variations of the flux density 
in HD 124314. They could be originated in changes in the 
mass loss rate, ionization or density structure that take 
place throughout the photosphere of the star. The corre- 
sponding time scale depends on the observing frequency. 
Scuderi et al. (1998) and Waldron et al. (1998) give ex- 
pressions to estimate the time scales of the variations. 

Let i?off be the characteristic radius of the emitting 
region, given by Wright & Barlow (1975) for a thermal 
outflow. When considering the mass loss rate predicted 
using Vink et al. (2000) models for HD 124314, R^s « 
1200 Rq and the transit time of the shock front in the 
radio regime will be > 5 days for variations in M . On the 
other hand, variations in ionization are characterized by 
time scales of > 20 days in the radio emission. The flux 
density changes must be faster towards higher frequencies, 
as i?cff diminishes. 

Following the model of Waldron et al. (1998), the time 
travel for a density disturbance to propagate from the X- 
ray radius to the radio radius is above 10 days. 

Taking the above considerations into account, is clear 
that the flux variation between CI and C2 can be ascribed 
to the mentioned phenomena. In contrast, it seems diffi- 
cult to find the same causes for the differences between the 
flux densities of the two observing runs during C2, which 
are separated just three days and occur at 6 cm but not at 
3 cm. A simpler explanation is to assume that the quoted 
flux density errors are underestimated. Nonetheless, it will 
be interesting to perform a series of observations with time 
scales of a few days towards this star, which is proposed as 
a counterpart to the unidentified gamma-ray source 3EG 
J1410-6147 (Romero et al. 1999). 



6.3. Comparison between radio-derived M and other 

estimates 

The mass loss rates derived here were compared to ob- 
served values found in the literature and from the new 
models of Vink et al. (2000). For the target stars, pre- 
vious estimates of M were derived from optical and UV 
profiles. The theoretical -expected- mass loss rates can be 
computed from Vink et al.'s models, as a function of the ef- 
fective temperature, luminosity, mass, and wind terminal 
velocity. For our calculations we have adopted the stellar 
parameters given by Vacca et al. (1996). One must bear in 
mind that predicted values were calculated with param- 
eters derived under a series of assumptions from models, 
depending on a certain spectral classification, which is not 
unique. 

Table 5 lists the mass loss rates obtained from the 

literature, Vink et al., and us. For the correlation we have 
assumed error bars of ~50% M in the values derived in 
this work. 

The result of the comparison showed no pattern. For 
the stars HD 94963, HD 97253, HD 101205, HD 135591 
and HDE 319718 the predicted and the radio values are in 
good agreement (see columns 3 and 4 of Table 5), although 
we would have expected radio signal around a 2a level for 
most of them. 

Despite HD 57060 spectral classification (Of+0), its 
predicted and optically derived mass loss rate, it was not 
detected at 3 or 6 cm. Using Eq. (1) we computed an ex- 
pected flux density, S'*, if M = Mopticai = 8.51 x 10~^ 
Mq yr-i. We found S^^^ « 0.6mJy > 4S'3cm and S^^^ « 
0.4mJy > 2S'6cm- The same relation of predicted and op- 
tically derived with radio mass loss rate happens to HD 
112244. Here S^^^ » 0.27mJy > 2S'.3cm. 

It is difficult to believe that the discrepancies between 
mass loss rates showed above might arise from underesti- 
mates in stellar distances, values which are, for both stars, 
widely accepted. Uncertainties in the wind terminal veloc- 
ities are evident from Table 1. Different methods of deter- 
mining yield to a range of values: ^ 1400 - 1800 km 
for HD 57060, and - 1600 2200 km s^^ for HD 112244. In 
our calculations we have adopted Woo(HD 57060) = 1800 
km s-i and ■Uoo(HD 112244) = 1900 km s'^ But the error 
on Voo does not suffice to account for the non-detections 
at radio continuum. We speculate that the presence of 
variability in the mass loss rates might be an explana- 
tion. Other possibility is that some of the assumptions 
undertaken in deriving M from optical lines (plus some 
theoretical modeling) (Hutchings 1976) are not entirely 
valid. 

For HD 135240, probably the expected mass loss rate 
computed using Vink et al. (2000) formulae is an overesti- 
mation. The radio upper limit is consistent with the mass 
loss rate from Hutchings (1976). 

The stars CD-47°4551, HD 124314 and HD 150136, 
with negative spectral indices, have theoretical mass loss 



Table 5. Non-radio (n-r), predicted (VKL), and radio- 
derived mass loss rates 




Star 


M (n-r) 


M (VKL)i 


M (radio) 




(Mo/yr) 


(Ms/yr) 


(Mo/yr) 


HD 57060 


8.51 10"" 


7.85 10"" 


< 3.03 10"" 




2.29 10"" 






CD-47° 4551 




1.81 10"^ 


2.39 10"^ 


HD 94963 




3.34 10"" 


< 7.16 10"" 


HD 97253 




5.22 10"" 


< 7.70 10"" 


HD 101205 




2.51 10"" 


< 8.38 10"" 


HD 112244 


5.01 10"" 


4.06 10"" 


< 2.87 10"" 


HD 124314 


7.08 10"" 


1.91 10"" 


1.81 10"^ 


HD 135240 


6.17 10"^ 


2.02 10"" 


< 2.03 10"" 


HD 135591 


3.98 10"^ 


1.81 10"" 


< 3.25 10"" 




5.50 10"^ 






HD 150135 




1.43 10"" 


5.89 10"" 


HD 150136 




5.07 10"" 


4.01 10"^ 


HD 163181 


7.08 10"" 


6.32 10"" 


2.37 10"" 


HDE 319718 




2.90 10"" 


< 4.62 10"" 



1: computed from Vink et al. (2000) 

rates lower than the values obtained from the radio ob- 
servations, as expected. The negative spectral indices are 
clear evidence of synchrotron contamination at 3 and 6 
cm, besides the free-free contribution. This contamina- 
tion may lead to the extreme case of an overcstimation 
in M of about an order of magnitude for HD 124314 and 
HD 150136, as Contreras et al. (1996) foimd for Cyg 0B2 
No. 9. The mass loss rate determined by Hutchings (1976) 
is similar to the radio one. 

In the case of HD 150135, {a > 0), the M derived 
here resulted greater than the predicted one, as for the 
stars with negative spectral index. In the model of White 
(1985) non-thermal emission is capable of producing a 
spectrum with a ~ 0.5. The search for variability can 
help to discriminate between a thermal or non-thermal 
emission mechanism. 

Bieging et al. (1989) reported that the observed mass 
loss rates were less than predicted for B supergiants of high 
luminosity. They pointed out that one source of error was 
that the model used to compute the expected rates ne- 
glected multiply scattered photons. Notwithstanding the 
model of Vink ct al. (2000) to derive M takes into account 
scattered photons, the radio-derived mass loss rate of HD 
163181 is less than predicted in a factor of 2. 

6.4- Origin of the non-thermal emission 

Regarding the non-thermal sources detected here, possi- 
ble origins of the emission are accretion of matter onto 
compact companions, shocks produced from line driven 
instabilities closer to the base of the wind, and in massive 
binary systems, at the wind collision region. Garmany et 



08 57 55.5 55.0 54.5 54.0 53.5 
RIGHT ASCENSION (J2000) 

Fig. 5. Contour map of CD-47° 4551 at 6 cm superposed 
on the J-band DSS image. 



al. (1980) have estimated that at most 6% of all O stars 
could have an (optically) undetected companion. Bieging 
et al. (1989) established that systems formed by an OB 
star and a compact companion comprises less than 1% 
of all luminous OB stars. Being the frequency of com- 
pact companions so low for early- type stars, the possibil- 
ity of detecting emission from this scenario becomes rather 
small. 

HD 150136 is a known binary star and HD 124314, a 
possible spectroscopic binary (see Table 1 for references). 
The determination of the spectral classification of the sec- 
ondary stars by means of optical observations is crucial 
to decide whether the non-thermal nature of these two 
candidates resides on colliding winds. 

CD-47° 4551 is not recognized either as a binary or 
multiple star but its radio image at 4.8 GHz exhibits a 
second source. The J-band image from DSS (see Fig. 5) 
shows that a second optical object could be present in the 
direction of this candidate. There exists the possibility 
that the weak radio source be originated by the presence 
of a companion of CD-47° 4551, separated ^ 9000 AU. 
A possible scenario could be the following: the strongest 
non-thermal source, positionally coincident with the star 
CD-47° 4551 could be generated by inner shocks in the 
wind of the Of star or shocks at wind collision region if 
the star has a close early-type companion. The second 
(and also non-thermal) source could be emission from a 
colliding winds region between the winds of CD-47°4551 
and a visual or wide-orbit object. 

The proposal stated by Chapman et al. (1999) and 
Dougherty & Williams (2000) that non-thermal emission 
could generally indicate wind interaction between a WR 
and a companion gives support to the idea of the existence 
of a windy companion for the present non-thermal sources, 
as Of stars are WR precursors. Optical studies of CD- 
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47° 4551 would be a fundamental tool to investigate its 
binary status. 

Unfortunately, the angular resolution of the current 
observations does not suffice to discriminate if the non- 
thermal emission comes from the wind of a star or from 
the site where two winds collide. This kind of investiga- 
tions must be carried out with very high angular resolu- 
tion telescopes, like MERLIN (see for example the results 
of Dougherty et al. 2000 related to WR 146) and they are 
limited by the source declination. In addition, variability 
studies would bring some information to help addressing 
this question, in the sense that if changes in the flux den- 
sity are present, the existence of an early type companion, 
and a wind-wind interaction region, is the most plausible 
explanation for the observed non-thermal radio continuum 
emission. 



7. Conclusions 

A number of conclusions can be drawn from the present 
observations: 

— Five out of thirteen OB stars were detected at 3 and 

6 cm wavelengths. The three brightest stars, CD -47° 
4551, HD 124314 and HD 150136, yielded negative 
spectral indices, giving < a > w —0.8, whereas HD 
150135 and HD 163181 gave positive spectral indices, 
< a > « +1.5. 

— None of the target stars, except HD 57060, had previ- 
ously been observed. This star presents now new (and 
lower) upper limits for the flux density (< 0.15 mJy 
at 3 cm) and for the mass loss rate (< 3 x 10~^ Mq 
yr""'^). A previous value of M, derived from optical 
data was 8.5x10"^ M© yr~^, and would correspond 
to a flux density of about 0.6 mJy at 3 cm. Variations 
(a decrease in this particular case) of the mass loss rate 
could explain the non detection above 0.15 mJy. 

— More than half of the detected stars show spectral 
indices consistent with contamination of non-thermal 
emission. If many surveys of continuum emission are 
taken into account, the results presented here increase 
the percentage of non-thermal stars up to 50% among 
detected 0-type stars. Non-thermal centimeter radio 
emission may be a common phenomenon in Of stars. 
A similar conclusion was reached by Chapman et al. 
(1999) concerning WR stars. These authors state that 
the emission would come in most cases from colliding 
winds. 

— Either shocks in the wind if the star has no early-type 
companion, or shocks in the region where winds collide 
can explain the non-thermal emission from three target 
sources. The stars could have a massive companion: 
HD 150136 is a binary star, HD 124314 is a suspected 
binary, and CD-47° 4551 presents both in radio and 
optical wavelengths a second source separated 5". 



— The flux density of HD 124314 was found to vary sub- 
stantially. Within two years the source almost dupli- 
cated its flux. 

— In this small sample, no clear correlation between de- 
tectability and distance, or spectral type, or luminosity 
class was found. However, the detected stars were the 
most luminous, early- type and closest candidates. 

— The theoretical predictions recently published by Vink 
et al. (2000) are considered a useful tool to estimate 
mass loss rates, and can also help to foresee in a general 
way which stars are more likely to be detected. Their 
predictions are in general agreement with the results 
presented here for undetected and thermal sources. 
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Table 1. Stellar parameters for the target stars 



Name 


Sp.Class. 


77) V 


(B V) 


d 

(kpc) 


Assoc. 




(km/s) 


Binary 
status 


HD 57060 


08.5If^ 
08Iaf+0'' 


4.90^ 


-0.14^ 


1.5 


NGC 23622 


-5.64^ 
-5.07^ 


1425*, 1800=5 
1700^ 


BIN'5 


CD-47° 4551 


05If^° 


8.39" 


0.89" 


1.7 






[1885] 




HD 94963 


06.5IH(f)i2 


7.18^ 


-0.09^ 


2.2 


Car 062^ 




[2545] 


C^ 


HD 97253 


05.5in(f)i2 


7.12^ 


0.15^ 


2.2 


Car OB22 




[2600] 
2740* 


bin2 


HD 101205 


07IIIn((f))i2 


6.48^ 


0.07^ 


2.5 


Cru OBl2 




SB2?2 


HD 112244 


08.5Iab(f)i2 


5.33^ 


0.01^ 


1.5 




-5.30® 


1950*', 1880^^ 
2160^,1575* 


SB1?2 


HD 124314 


06V(n)((f))i^ 


6.64^ 


0.21^ 


1.0 




-5.15® 


2500* 


SB1?2 


HD 135240 


07.5in((f))i5 


5.08^ 


0.06^ 


1.0 




-6.21^" 


2460*, 2700=5 


bin2 


HD 135591 


07.5IH((f))i5 


5.43^ 


-0.09^ 


1.2 




-6.40^-7.26^'^ 


2600^'^, 2235* 


BIN^® 


HD 150135 


06.5V((f))i5 


6.89^ 


0.13^ 


1.4 


Ara OBla^ 




[2455] 
3160* 


SB1?2'6 


HD 150136 


05mn(f)is 


5.64^ 


0.14^ 


1.4 


Ara OBla^ 




BIN« 


HD 163181 


Bllape^" 
BN0.5Iap2i 


6.43i« 


O.52I8 


1.4 




-5.15® 


520* 


BIN^" 


HDE 319718 


07II122. 03111^=' 


10.4322 


1.4822 


1.7 


NGC 63572-^-24 




[2295] 





1 - Cruz-Gonzalcz ct al. (1974); 2 Gics (1987); 3 Drochscl ct al. (1980); 4 Howartli ct al. (1997); 

5 - Hutchings & von Rudloff (1980); 6 - Chlcbowski & Garmany (1991); 7 Garniany, Couti & Masscy (1980); 

8 - Hutchings (1977); 9 - Hutchings (1976); 10 - Garrison et al. (1977); 11 - Drilhng (1991); 12 Walborn (1973); 

13 - Vijapurkar & Drilling (1993); 14 - Bernabeu et al. (1989); 15 - Walborn (1972a); 16 - Prinja & Howarth (1986); 

17 - Conti & Garmany (1980); 18 - Lindroos (1985); 19 - Garmany, Conti & Chiosi (1982); 20 - Levato et al. (1988); 

21 - Walborn (1972b); 22 - Neckel (1984); 23 - Johnson (1973); 24 - Felli et al. (1990) 

Table 2a. Positions and flux densities of the detected stars 



Star Optical positions Radio positions tint Sscm Seem 
R.A.(2000) Dec.(2000) R.A.(2000) Dec.(2000) (min) (mjy) (mjy) 



CD-47° 4551 


08 


57 


54.62 


-47 44 


15.7 


08 


57 


54.61 


-47 


44 


15.7 


160 


1.77±0.05 


2.98±0.05 


HD 124314 


14 


15 


01.61 


-61 42 


24.4 


14 


15 


01.61 


-61 


42 


24.2 








CI: 02/16/1998 
























150 


1.88±0.05 


2.72±0.08 


C2: 03/30/2000 
























195 


2.71±0.08 


4.50±0.07 


C2: 04/02/2000 
























130 


2.88±0.05 


3.78±0.05 


HD 150135 


16 


41 


19.45 


-48 45 


47.6 


16 


41 


19.45 


-48 


45 


47.8 


185 


0.28±0.03 


0.10±0.04 


HD 150136 


16 


41 


20.41 


-48 45 


46.6 


16 


41 


20.42 


-48 


45 


46.8 


185 


2.61±0.03 


5.57±0.03 


HD 163181 


17 


56 


16.08 


-32 28 


30.0 


17 


56 


16.07 


-32 


28 


29.5 


135 


0.44±0.05 


0.20±0.05 



Table 2b. Flux density limits of the undetected stars 



Star 


Optical 


positions 


^int 


'S'scm Seizin 


R. A. (2000) 


Dec. (2000) 


(min) 


(mjy) (mjy) 


HD 57060 07 18 40.38 


24 33 31.3 


180 


<0.15 <0.15 


HD 94963 10 56 35.79 


61 42 32.4 


140 


<0.14 <0.17 


HD 97253 11 


10 42.02 


-60 23 04.3 


140 


<0.15 <0.15 


HD 101205 11 38 20.36 


-63 22 21.9 


145 


<0.13 <0.20 


HD 112244 15 55 57.13 


-56 50 08.9 


140 


<0.13 <0.15 


HD 135240 15 16 56.90 


60 57 26.1 


170 


<0.13 <0.16 


HD 135591 15 18 49.14 


60 29 46.8 


180 


<0.15 <0.15 


HDE 319718 17 24 42.90 


-34 11 48.0 


180 


<0.15 <0.40 


Table 4. Ra( 


lio-(l(;rivo(l 


mass 


loss rates 






Star 






Z 7 


53cm 


log(M) 




(K) 








(Mo/yr) 


HD 57060 


14300 


1.5 


1.0 1.0 


5.1 


< -5.52 


CD 47° 4551 


17900 


1.5 


1.0 1.0 


5.3 


(-4.62) 


HD 94963 


16500 


1.5 


1.0 1.0 


5.2 


< -5.15 


HD 97253 


17600 


1.5 


1.0 1.0 


5.2 


< -5.11 


HD 101205 


16000 


1.5 


1.0 1.0 


5.2 


< -5.08 


HD 112244 


13700 


1.5 


1.0 1.0 


5.0 


< -5.54 


HD 124314 


17400 


1.3 


1.0 1.0 


5.2 


(-4.74) 


HD 135240 


15400 


1.5 


1.0 1.0 


5.1 


< -5.69 


HD 135591 


15400 


1.5 


1.0 1.0 


5.1 


< -5.49 


HD 150135 


16900 


1.5 


1.0 1.0 


5.2 


-5.23 


HD 150136 


18200 


1.5 


1.0 1.0 


5.3 


(-4.40) 


HD 163181 


11600 


1.6 


0.9 0.8 


5.0 


-5.63 


HDE 319718 


15900 


1.5 


1.0 1.0 


5.2 


< -5.34 



